Chapter 1 

The initial conditions of star formation: cosmic 
rays as the fundamental regulators 

Padelis P. Papadopoulos and Wing-Fai Thi 



Abstract Cosmic rays (CRs) control the thermal, ionization and chemical state of 
the dense H2 gas regions that otherwise remain shielded from far-UV and opti- 
cal stellar radiation propagating through the dusty ISM of galaxies. It is in such 
CR-dominated regions (CRDRs) rather than Photon-dominated regions (PDRs) of 
H2 clouds where the star formation initial conditions are set, making CRs the ulti- 
mate star-formation feedback factor in galaxies, able to operate even in their most 
deeply dust-enshrouded environments. CR-controlled star formation initial condi- 
tions naturally set the stage for a near-invariant stellar Initial Mass Function (IMF) 
in galaxies as long as their average CR energy density Ucr permeating their molec- 
ular ISM remains within a factor of ~10 of its Galactic value. Nevertheless, in the 
extreme environments of the compact starbursts found in merging galaxies, where 
UcR^(few) x 10 3 UcR.Gai, CRs dramatically alter the initial conditions of star forma- 
tion. In the resulting extreme CRDRs H2 cloud fragmentation will produce far fewer 
low mass (<8 M E ) stars, yielding a top-heavy stellar IMF. This will be a generic fea- 
ture of CR-controlled star-formation initial conditions, lending a physical base for 
a bimodal IMF during galaxy formation, with a top-heavy one for compact merger- 
induced starbursts, and an ordinary IMF preserved for star formation in isolated 
gas-rich disks. In this scheme the integrated galactic IMFs (IGIMF) are expected 
to be strong functions of the star formation history of galaxies. Finally the large, 
CR-induced, ionization fractions expected for (far-UV)-shielded H2 gas in the CR- 
DRs of compact starbursts will lengthen the ambipolar diffusion (AD) timescales 
so much as to render the alternative AD-regulated rather (Jeans mass)-driven star 
formation scenario as utterly unrealistic for the ISM in such galaxies. 
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1.1 Cosmic rays versus far-UV and optical photons 

Much of the interstellar medium (ISM) in galaxies resides in the so-called Photon- 
dominated regions (PDRs) where stellar far-UV photons (6eV</jV<13.6eV) de- 
termine the thermal, ionization and chemical state of atomic (HI) and molecular 
(H2) hydrogen gas |T). For HI these phases are the Warm Neutral Medium (WNM, 
with T kin ~10 4 K and n~(0.1-l)cirr 3 ) and the Cold Neutral Medium (CNM, with 
T kin ~(80-200)K and n~(10-100)cm- 3 ) [2]. For the molecular PDRs in H 2 clouds 
on the other hand: n~(10 2 -10 6 )cirr 3 and T kin ~(30-100)K. Such PDRs start from 
the surface of molecular clouds (MCs) with a layered distribution (Fig. Il.lt of 
atomic and ionized species on the outside (HI, C + , C, O), and continue deep in- 
wards where molecules such as CO dominate EH), marking the location of such 
clouds via their mm/submm rotational line emission. The thermal, ionization, and 
chemical state of all these phases is determined by the far-UV photons whose influ- 
ence, especially on the ISM chemistry, can extend very deep inside the low-density 
MCs (n(H 2 )~(100-500)cirr 3 ) typical in the Galaxy. The dust in PDRs is heated 
via absorption of the ambient far-UV and optical radiation and cools via continuum 
thermal IR/submm emission while the gas is heated via electrons ejected by the 
photolectric effect on the concomitant dust grains and polycyclic aromatic hydro- 
carbons (PAHs). Despite the gas receiving ~10 2 -10 3 less heating per unit volume 
than the dust, in PDRs it is always T k i n >T ( j ust . This is because of the much less ef- 
ficient cooling of the gas via spectral lines while dust cools via continuum radiation 
across a wide range of wavelengths (typically peaking at IR/submm wavelengths). 
Turbulent heating of the gas in MCs will further reinforce the T^i^Tdust inequality, 
since it does not involve the dust reservoir. 



Photodissociation Region 



Fig. 1.1 The structure of a 
typical PDR/CRDR. Atomic 
and ionized species reside 
on the (far-UV)-illuminated 
outer parts (left) while neutral 
atoms (O) and molecules 
such as CO dominate deeper 
inside, the domain of CRDRs 
and CR-controled thermal and 
chemical gas states. 




For MCs in the typical ISM conditions of spiral disks where (A v )^3-7 (in units 
of optical extinction), there can only be very little H 2 gas residing in their dense (far- 
UV)/optical-shielded regions. Its mass fraction per MC is ~(l-5)%, a census found 
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from surveys of CO J=l-0 (total H2 mass indicator) and rotational transitions of 
dense gas tracer molecules such as the HCN (mostly J= 1—0), and using the HCN/CO 
line ratio to determine fdense = 

M(n > 10 4 cm- 3 )/M tot (H 2 ) Eg). Here is important 
to note that the more complex heavy rotor molecules such as HCN will necessarily 
trace the dense but also (mainly) the (far-UV)-shielded mass of MCs since, unlike 
H2 and its tracer molecule CO, they have much lower dissociation potentials (—5 eV 
for HCN versus -13 eV for H 2 , and -lOeV for CO). 

Simple approximate estimates of the mass fraction contained in the PDRs of (far- 
UV)-irradiated MCs can be found by computing the mass of outer atomic HI layer 
marking the HI— >H 2 transition, itself comparable to the mass of the warm PDR H 2 
gas layer extending further inwards. Thus the PDR-residing gas column density (HI 
and H 2 ) will be — 2xN tr (HI). The latter can be computed for a given radiation field 
Go (in Habing units), metallicity Z (Z=l being solar), average gas density n, and an 
H 2 formation rate Rf. In units of optical extinction this transition layer is ([9 1): 



A? = 1.086§^yln 



1 + ^r7* 



(i.i) 



where k D =4x 10~ n s -1 is the H 2 dissociation rate (forG =l), Rf— 3x 10~ 17 cm~ 3 s -1 
its formation rate on grains (for typical CNM HI), and 4>=6.6x 10 _6 v / ^Z 1 / 2 ^ FUV 
is the H 2 self-shielding function integrated over the HI/H 2 transition layer, with 
<sFUV=Ofuv / CTy— 2-3 being the dust cross section ratio for far-UV and optical light. 
Using existing formalism (|9|) the PDR-related gas mass fraction per MC then is: 



fpDR ~2X 




(1.2) 



assuming spherical and uniform MCs that do not cross-shield and hence each re- 
ceives the full (far-UV)/optical radiation field, assumptions that make the computed 
fpDR a maximum. The average optical extinction (A v ) per MC is of course finite and 
a function of the ambient ISM conditions, mainly the gas pressure (dominated by su- 
personic non-thermal rather than thermal velocity fields). Expressing an empirical 
Galactic (average MC density)-(size) relation: n«:(2R) , along with its expected 
normalization in terms of cloud boundary pressure P e (|9 |) in terms of average op- 
tical extinction yields: 

< A ->~°- 22z =Gra)" 2 - (i - 3 > 

where n — 1500cm -3 . For average MC densities of —100cm -3 , Solar metalic- 
ities, ^puv=2, and G o =5-10 (for typical star-forming gas) Equation 1.1 yields 

Av tt ^~0.50-0.75, while for typical interstellar pressures in galactic disks P e /kB— (1- 
2)xl0 4 Kcm -3 it would be (A v )-3.3-4.7. The latter is the A v range over which 
the abundance of C + plummets, giving rise to C and CO while the free electron 
abundance is much reduced (see Fig. II. 2b . 
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Fig. 1.2 The PDR— >CRDR transition zone as marked by the rapidly decreasing C+ and e~ abun- 
dances deep inside a (far-UV)-illuminated molecular cloud (radiation field incoming from the left). 
Note the dramatic fall of the electron abundance beyond A v ~4. 



For the finite MCs this outer HI/(warm H 2 ) PDR gas still contains f PDR — 0.74-1 
of their mass (from Equation 1.2), i.e. most of the HI and H2 gas in typical star- 
forming ISM environments resides in PDRs. Later we will see that this will not be 
so in the high-pressure ISM of Ultra Luminous Infrared galaxies (ULIRGs), merger 
systems with very IR-luminous (Lir>1O 12 L ) compact starbursts (Sect. [L4"l ). 



1.2 The physical conditions of the gas: PDRs versus CRDRs 

The wide range of physical conditions of the HI and H2 gas in PDRs (Sect. 11.11 ) is 
the direct result of: a) the widely varying strengths of the dust-absorbed far-UV and 
optical radiation (and thus of the gas photoelectric heating per H, and dust heating) 
within the PDRs, and b) the large range of densities (and thus of the cooling power 
per H) that are expected in the supersonic ally turbulent MCs (|7|). Dissipation of 
the supersonic turbulence via shocks (J8)) will further contribute to the wide range 
of temperatures observed in the PDR-residing gas of MCs. Moreover their average 
electron abundance x(e)=n e /[2n(H2)] changes by ^2-3 orders of magnitude, start- 
ing from a maximum of ^3 x 10~ 4 (that corresponds to a fully ionized carbon). 



1 The initial conditions of star formation: cosmic rays as the fundamental regulators 



5 



The dense gas cores (n~(10 5 -10 6 )cm~ 3 ) in Galactic CRDRs on the other hand 
span a much narrower range of temperatures and ionization states, with (Ty n )^10 K 
(and ^30% dispersion), and x(e)~(few)x 10~ 8 (varying by less than a factor of 
~5). These CRDR conditions, widely observed in the Galaxy ( iflOl . JTT|), are a 
result of a nearly complete lack of far-UV and optical photons in such regions (thus 
no photolectric gas heating or strong dust heating), nearly complete dissipation of 
supersonic gas motions ([12]), and the onset of strong gas-dust thermal coupling 
setting Tki n ~Tdust- This simply means that cold dust grains inside CRDRs, warmed 
only by the feeble IR radiation able to leak inside them, can now act as powerful 
'thermostats' of the H2 gas. The latter is now heated almost exclusively by primary 
and secondary ionizations induced by low-energy CRs (10MeV<EcR<100MeV) 
|[T3l which, unlike far-UV and optical photons, penetrate deep inside such regions, 
and thus fully regulate their thermal, ionization and chemical states. 

The aforementioned CR-controlled dense H2 gas states in ordinary star-forming 
environments is expected also on broad theoretical grounds, namely: a) high density 
gas regions will be typically nested well inside less dense gas structures of super- 
sonic MCs [14], and thus will always 'see' a much-attenuated stellar radiation field, 

b) high-density gas is also strongly cooling gas (since (line-cooling)^ [n(H2)] 2 ), and 

c) turbulent energy injected on the largest gas structures with low average densities 
~(10-10 3 ) cm~ 3 will have fully dissipated in the high density gas of CRDRs. It is 
actually these three reasons, along with the strong thermal coupling of the gas in 
CRDRs to a necessarily cold dust (for lack of far-UV/optical photons) reservoir, 
that make it very difficult for large masses of warm and dense H2 gas to exist in 
ordinary star-forming galaxies. Finally, as already mentioned in 11.11 only few% of 
the total H2 gas mass resides at such high densities anyway, as expected also from 
MC simulations ([7]) under conditions found in ordinary star-forming spiral disks. 



1.2.1 Computing T^[ n andx(e) in CRDRs: a simpler method 

In CRDRs the Tj™ n) and x min (e) minimum values expected for H2 gas in galaxies 
are set. However determining them involves solving the coupled equations of ther- 
mal, chemical and ionization balance in such regions since CRs heat the gas but also 
affect the abundance of various molecular and atomic coolants. We can nevertheless 
make some simplifying assumptions that will keep the physics involved transparent 
while demonstrating the range of Tj.™ and x m i n (e). Following standard methods 
(El) T^ n ' in CRDRs can be computed by solving the equation of thermal balance: 

r C R = Aline +A g _ d , (1.4) 

where rcR oc CcR n (H2) is the CR heating with ^cr^ Ucr being the CR ionization 
rate, and Ucr the average CR energy density. For Tki n <50K (which encompasses 
what is expected for CRDRs of ordinary star-forming galaxies), the line cooling 
term Ai; ne is dominated by CO lines. However for the much more extreme CRDRs 
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expected in ULIRGs, higher temperatures are possible and thus cooling from the 
atomic fine structure lines OI (at 63 jim) and C + (at 158 /im) can become impor- 
tant. The C + abundance in CRDRs is not necessarily negligible as it is no longer 
controlled by a stellar far-UV radiation field (as in PDRs), but by ISM-CR interac- 
tions inducing internal UV radiation which destroys CO and produces C + , a clear 
demonstration of the coupled thermal and chemical states of the gas in CRDRs. Fi- 
nally the term A g _d denotes the all-important (in CRDRs) cooling of the gas via 
the gas-dust coupling which depends on the gas and dust temperatures. Thus the 
thermal balance equation becomes 

Aco(Tkin) + A g _d(Tkin,Tdust) + Aoi63(Tkin) +A c +(Tki n ) = Tcr, (1.5) 

(see the Appendix for the detailed expressions used). We have not considered tur- 
bulent gas heating since we want to compute the minimum CR-controlled gas tem- 
peratures and since supersonic turbulence (and thus shock-heating) has typically 
dissipated in the dense pre-stellar gas regions inside CRDRs. 




1 I ^ ^ ^ I 

10 3 10 4 10 5 10 6 

Gas density (cm 3 ) 

Fig. 1.3 The equilibrium temperature of (far-UV)-shielded regions in CRDRs computed from 
equation 1.5. The CR energy densities range from those expected in ordinary star-forming en- 
vironments Galaxy [Ucr=(0.2-10)xUcr Gal] to those expected in the extreme CRDRs of ULIRGs 
[UcR=(10 2 -10 4 )xUcR.Gal]- A dust temperature of T ( j ust =8 K, and a Galactic CR ionization rate of 
Ccr=5x 10~ 17 s have been assumed. 



1 The initial conditions of star formation: cosmic rays as the fundamental regulators 



7 



In order to solve the gas thermal balance equation in principle we also must solve 

for Tdust in CRDRs. For the purposes of computing T^ 1 ' we set Tdust=8K, which 
is typical for CR-dominated cores in the Galaxy. Any stronger IR radiation field 
leaking inside such cores (as expected in starbursts) or remnant turbulent gas heating 
can only raise the gas temperatures computed here. In Fig. 11.31 the T^™"' values 
show that CR-heated dense gas (>10 4 cm~ 3 ) remains cold (~(10-15)K) for Ucr 
within ^(0.2-10)xUcR.Gai, rising decisively only when UcR~(few)x 10 3 UcR,Gai 
where the thermal balance equation yields Tj,™ n '~(40-100)K. 

Solving for the coupled chemical and thermal states of dense gas in CRDRs 
does not change much of the aforementioned picture while it makes clear another 
important aspect that of CR-controlled thermal states of the gas, namely that high 
CR energy densities will induce a strong gas-dust thermal decoupling, with 
remaining significantly higher than T(j ust , even at high gas densities (Fig. |1.41 >. 




Fig. 1.4 The gas and dust temperature profiles (solid and dotted lines of similar color) versus opti- 
cal extinction inside a molecular cloud. These were computed by solving the coupled thermal and 
chemical balance equations (|19|), while also assuming a far-UV radiation field (incident on the 
A v =0 surface of plane clouds) that scales in an identical fashion as the average CR energy density 
pervading the H2 cloud. Their common scaling runs as (l,10,10 2 ,10 3 ,10 4 )xGalactic. Unlike in 
Fig. l 1 ,3l the dust is significantly wanned by the stronger (far-UV)/optical radiation fields absorbed 
by the outer cloud layers and the re-radiated dust IR emission that penetrates deeper into the cloud. 
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This happens even as the dust temperatures deep inside H2 clouds rise because 
of an incident radiation that is now set to scale identically to the average CR energy 
density (both assumed proportional to the average star-formation rate density p s fr)- 
This CR-induced thermal decoupling persists deep inside the cloud (A v =5-10) even 
at the highest gas densities explored here where gas-dust thermal interaction will 
be the strongest. It must also be noted that for a given average p s f r in a galaxy, 
the resulting gas-dust temperature inequality will likely be even larger since cross- 
shielding between H2 clouds as well as dusty WNM and CNM HI gas interdispersed 
between such clouds will reduce the average radiation field incident on them, but not 
the deeply penetrating CRs. Thus while Ucr k Pj/>, the scaling of average G with 
p s f r will be weaker than linear (unlike what was assumed in Fig. ll.4l for simplicity). 

From Fig. |1.4| it can also be discerned that even for strongly boosted far-UV radi- 
ation fields, gas temperatures are solely CR-controlled and remain almost invariant 
with A v beyond Ay ~5. This characteristic of CRDRs is of great importance when 
it comes to the initial conditions of star formation, and the near-invariant mass scale 
of the stellar IMF found in the Galaxy and ordinary star-forming galaxies. 

An order of magnitude estimate of the CR-regulated equilibrium ionization frac- 
tion of (far-UV)-shielded gas can also be found without resorting to solving the full 
chemistry network involved. Following an analytical approach (|20|) yields 



x(e) = 2x 10" 'r, 



-7-1 ( "ch \ 



1/2 



\ 1/2 / \ 1/2 

8n(H 2 ); V8n(H 2 ); 



(1.6) 



where n c h^500 yr^ d 17 J cm is a characteristic density encapsulating the effects 
of cosmic ray ionization and ambient metallicity on the ionization balance (r^ is 
the normalized gas/dust ratio r g( j=[(G/D)/100] with G/D(gas-to-dust mass)=100 for 
Solar metallicities. The main assumptions underlying the validity of Equation 1 .6 
are three, namely: a) all the gas is H2 and all the carbon is locked in CO (so C + does 
not dominate the ionization), b) all negative charge is carried by the free electrons 
and molecular ions are destroyed mostly by recombination with these free electrons, 
and c) the metal abundance is small with respect to the heavy molecule abundance. 

For the Galaxy (£cR.Gai=5x 10" 17 s 1 ) it is n c h=2.5x 10 3 cm" 3 , and for a typical 
dense core of n(H2)=10 5 cm" 3 it would be x(e)Gai~2.4x 10" 8 . A boost of £cr by a 
factor of 10 yields x(e)=8.4x 10" 8 , retaining the low ionization fraction necessary 
for an ambipolar-diffusion regulated star-formation to proceed in such cores ( [ 20 1). 



1.3 CRDRs and the star formation initial conditions in galaxies 

CRDRs are the ISM regions where the initial conditions of star-formation are set. 
Thus their thermal and/or ionization state is expected to play a critical role on the 
emergence of the stellar IMF and its mass scale. This remains the case for both main 
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star formation scenaria proposed, namely: a) star formation driven by gravitational 
collapse of high overdensity peaks in supersonically turbulent H2 clouds, and b) 
regulated by ambipolar-diffusion (AD) in dense cores whose low x(e) allows the 
magnetic field lines to slip away making them supercritical (i.e. M core /M<j>> 1 where 
M<j,=0.12<I>/G 1 ' /2 , and <P is the magnetic field flux threading the cores, [21]) and 
able to collapse towards star formation. The only difference will be which CRDR 
gas property: the Jeans mass or x(e), is important in determining the stellar IMF. 

It is worth emphasizing that star formation initial conditions set in CRDRs natu- 
rally explain the first mystery regarding the stellar IMF, namely its near-universality 
in ordinary ISM environments. Indeed the large range of thermal, turbulent, and 
ionization states present for the H2 gas phase in the Milky Way, makes the emer- 
gence of a near-universal stellar IMF an unlikely outcome even for the Galaxy if 
such states were all used as star formation initial conditions. The gas in CRDRs 
on the other hand, lying beyond the reach of varying far-UV radiation fields and 
turbulence-induced shock heating, with a near-isothermal state at TW^T^st, near- 
thermal motions, and the lowest possible x(e) values, is a natural choice of well- 
protected initial conditions of star formation that can then lead to a near-invariant 
stellar IMF mass scale. Similar arguments have been made in the past ([22]), but 
only recently it was demonstrated just how robust the cold thermal state of CRDRs 
remains irrespective of the average conditions prevailing at the outer boundaries of 
the MCs that contain them (OH). Finally, as long as the average CR energy density 
pervading CRDRs remains within a factor of ^10 of its Galactic value, neither the 
temperature nor the ionization fraction within CRDRs changes substantially. This 
state of affairs is however dramatically altered in the compact, merger-induced star- 
bursts found in ULIRGs. 



1.4 CRDRs in compact starbursts: new initial conditions for star 
formation 

Some of the most spectacular star formation events in the local and the distant Uni- 
verse are found in strong galaxy mergers. In such systems tidal torques funnel the 
bulk of the molecular gas of their gas-rich disk progenitors into regions not much 
larger than the Galactic Center (EH). There they settle into small (2r<(100-300)pc) 
and very turbulent H2 gas disks, fueling extreme star-forming events ([25, 24]). 
The H2 gas in ULIRGs, with one-dimensional Mach numbers of MulirGs~(5- 
30) x M S pi ra i s , is expected to have larger mass fractions at densities of >10 4 cm -3 
(0), and this is indeed observed (|5]|6]). 

However most of the large dense H2 masses found in ULIRGs cannot reside 
in PDRs, a result of the high (turbulent) pressures (P e /kB~10 7 )Kcm~ 3 ) and high 
gas densities in these galaxies. Indeed, from Equation 1.3 (and a Solar metallicity) 
such high gas pressures yield (A v )^100, while for the high average gas densities 

of >10 4 cm- 3 typical for ULIRGs (Eg)), it is Ai tr) <2 (Eq. 1.1) even for far-UV 
radiation fields as strong as G o =10 4 . Thus any warm gas mass fraction residing in 
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PDRs will be < 15% (Equation 1 .2). In summary, while much larger dense gas mass 
fractions are expected and indeed observed in ULIRGs than in isolated star-forming 
gas-rich spirals, most of their dense gas mass will not reside in PDRs. Neverthe- 
less recent observational results indicate large amounts of warm and dense gas in 
some ULIRGs, with average gas temperatures much larger than those deduced for 
the concomitant dust ([27]). Since a massive warm and dense H2 phase cannot be in 
PDRs, far-UV/optical photons are inadequate energy suppliers of its observed state 
in ULIRGs. Only high CR energy densities (and/or strong turbulent heating) is ca- 
pable of volumetrically heating the H2 gas (unlike far-UV photons in PDRs) while 
also maintaining a significant Tki n >Td U st inequality. 

CR energy densities of UcR^[(few)x 10 2 -10 4 ] xUcR.Gai are expected in the 
compact starbursts found in ULIRGs ([29 1) and, as indicated by Fig. 11.41 such high 
Ucr can easily maintain high gas temperatures while the concomitant dust remains 
much cooler. Such conditions have been inferred in the past for the Galactic Cen- 
ter ( POl ) as well as for the center of NGC253 (GUIS)). However in the case of 
ULIRGs these conditions prevail over galaxy-size H2 gas reservoirs that are ~2-3 
orders of magnitude more massive. Thus temperatures of Tki n ^(50-100)K are the 
new minimum values possible for gas in ULIRGs, and represent dramatically dif- 
ferent initial conditions for star formation in their ISM environments. 

The thermal state of dense gas in CRDRs is thus no longer a near-invariant but 
strongly altered by the high CR energy backgrounds expected in ULIRGs. In the so- 
called gravoturbulent H2 cloud fragmentation scenario ([33 1) this large temperature 
boost across a wide density range (Fig. 11.31 is expected to surpress low mass star 
formation. This is because the Jeans mass of the CR-heated gas in CRDRs, 



M 



(c) _ / k B T k 



Jeans 



\GjUm H2 



3/2 



^3/2 


' n c (H 2 ) " 




10 4 cm- 3 _ 



n -1/2 



Ms, (1.7) 



changes from M^~0.3Mg for n c (H 2 )=10 5 cirr 3 , and T k =10K (typical for the 

Galaxy) to M^ ns ~(3-1O)M for the same density and ^=(50-1 10) K expected 
to be the new minimum in the extreme CRDRs of ULIRGs. It is worth pointing 

(c) 

out that an Mj e ^ ns ^0.3M e nicely matches the characteristics mass of the stellar 
IMF (the so-called IMF 'knee')- Its observed near invariance under a wide range 
of ISM conditions (e.g. radiation fields) has been demonstrated theoretically (|34|), 
leaving CR energy densities as the sole star-formation feedback mechanism that can 
decisively alter the initial conditions of star formation in galaxies. This mechanism, 
driven by the massive stars (whose strong winds and explosive deaths provide the 
CR acceleration), will operate unhindered in the dusty ISM environments of extreme 
merger/starbursts (unlike photons), and it can reach where supersonic turbulence 
cannot, the small, dense and (far-UV)-shielded cores deep inside H2 clouds. 

In Fig. 1 1 .51 the Jeans masses are shown for a density range encompassing that 
of typical star-forming H2 clouds, demonstrating that for UcR/UcR,Gai>10 3 , Mj ea n S 
increases by a factor of ^10 across the entire density range. This will invariably 
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Gas density (cm' 3 ) 

Fig. 1.5 The Jeans mass estimated from Equation 1.7 and the gas temperatures inside CRDRs, 
which were determined by solving the fully coupled thermal and chemical balance equations. The 
temperatures shown are averages of the A v =5-10 part of the H2 clouds where CR become dominant 
regulators of the gas thermal and chemical states (Fig. 11.4b . 

lead to higher characteristic mass M^. h for newly formed stars, and thus a top-heavy 

stellar IMF (J22] [34)). For an M^~M jeans (n c ) at the typical density of Galactic 
CRDR dense gas cores (n c ~(few)x 10 s cm~ 3 ) Fig. 11.51 also shows that variations 
of UcR/UcR.Gai-0.2-10 would leave M^ ) ~(O.3-O.5)M and still compatible with 

a standard stellar IMF. This is important since a larger sensitivity of M^. h on Ucr 
would have been incompatible with its observed near-invariance in ordinary star- 
forming environments (where variations of the star-formation rate density p s f r , and 
thus of Ucr, by factors of ^5-10 are not extraordinary). 

It must be noted that the choice of the particular (n,Tk) phase whose Mj eans 
sets M^. h is not without ambiguities even if some concensus has been established 
on the decisive role of Mj eans (r,t) in driving the outcome of the gravoturbulent 
process at each spatial and temporal point (r,t). Arguments have been made both 
for M^°c(Mjeans) at the onset of a gravoturbulently driven MC fragmentation pro- 
cess ( IT3311 ). and for M^~M jeans ( n c,Tk) at some characteristic gas density n c . For 
gravoturbulent-regulated star-formation n c is the density where efficient gas-dust 
thermal coupling of (far-UV)/optically-shielded gas lowers its temperature to the 
minimum value possible, i.e. that of ~Td us t, rendering such gas regions nearly 
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isothermal. For the Galaxy this occurs at n c ^10 5 cm~ 3 , which is also similar to 
the density of regions where the transition from supersonically turbulent to near 
thermal gas motions occurs ( |[T0ll22l '). 

For CRDRs with high average Ucr like those expected in compact starbursts, 
Fig. 1 1.41 makes obvious that gas-dust coupling is no longer adequate to lower T kin 
to that of the concomitant dust for densities as high as 10 6 cm~ 3 . Thus a charac- 
teristic gas density can no longer be defined by the onset of thermal equilibrium 
between gas and dust in CRDRs, at least within the range of densities where ap- 
preciable amounts of mass can be found in ordinary MCs. This leaves the choice 
of M^—M jeans (n c ,T]c) open, and thus it is worth examining the possible direc- 
tion its value would take in extreme CRDRs. In the CR-innudated H2 clouds in- 
side compact starbursts n c can now be set only as the typical density of regions 
where supersonic turbulence has fully dissipated and near-thermal motions domi- 
nate. The average density of such transition regions can be determined from two 
well-known scaling relations found for H2 clouds: OV(r) =(7 °( r /P c ) h (linewidth-size) 
and (n)=n (r/pc) _1 (density-size) after setting (7v(min)=(3kBTi c //xmH 7 ) 1,/2 as the 
minimum linewidth possible, and solving for the corresponding mean density, 

(nc)=n °bk^rJ ~ lln °(km^) UkJ ■ (L8) 

For h=l/2 (expected for pressurized H2 virialized clouds l35ll ). <7 =1.2kms _1 and 
n =(few)x 10 3 cm~ 3 (from observations), it is (n c )~(few)x 10 4 cm~ 3 . For the CR- 
boosted T kin (min) of extreme CRDRs T k (min)~(50-100)K: (n c )~fewx 10 3 cirr 3 . 
Thus the (turbulent gas)— >(thermal core) transition in extreme CRDRs will occur at 
lower densities as well as higher (CR-boosted) values than in ordinary Galac- 
tic ones. This will then shift a M^~Mj ean s(nc>Tk) both vertically (higher T k i n ), as 
well as to the left (lower n(H2j) of Fig. 11.51 yielding even higher M^. h values than 
those resulting from only a vertical shift of M] ean s because of CR-boosted gas tem- 
peratures. Thus irrespective of the H2 cloud fragmentation details, and as long as 
the Jeans mass plays a decisive role in it, a much larger characteristic mass M J is 
expected for stars forming in the CRDRs of extreme starbursts. 



1.4.1 Ambipolar-diffusion regulated star formation 

The CR-induced boost of the ionization fraction x(e) inside the CRDRs of extreme 
starbursts dramatically alters the star-formation initial conditions relevant for the 
ambipolar-diffusion (AD) regulated star formation scenario. Indeed, the much larger 
ionization fractions that can now be achieved deep inside H2 clouds can in prin- 
ciple: a) keep the magnetic field lines strongly "threaded" into molecular gas at 
much higher densities, and b) as a result render much of its mass incapable of star- 
formation, at least in the simple (AD)-regulated star formation scenario. This stems 
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from the now much longer ambipolar diffusion timescales: Tad~1-6x 10 14 x(e)yrs 
needed for a CR-ionized core with density n(H2) to lose magnetic flux and collapse 
(|20|). Indeed, after inserting x(e) (from Eq. 1.6) it is: 



- = 3.2xl0\ d (^) 



1/2 



\ 1/2 / \ 1/2 

nch\ / n ch \ 
+ 8n(H 2 ); \8n(H 2 )J 



yrs. (1.9) 



In Fig. 11.61 Tad is shown as a function of the CR ionization rate £cr (assumed 
^Ucr) for three typical densities of cloud cores inside CRDRs. It is obvious that the 
much higher £cr values expected in the CRDRs of (U)LIRGs will make the corre- 
sponding Tad timescales much longer than those expected in the Galaxy. Moreover 
Tad m (U)LIRGs will be longer than their so-called gas consumption timescales 
f C ons = M(H2)/SFR, which are observationally determined using M(H2) (usually ob- 
tained via CO J=l— 0) and the star formation rate (SFR) (usually obtained from Ljr). 

For any theory of star formation T cons represents a natural maximum for any the- 
oretically deduced characteristic star-formation timescale T* c h since star-formation 
feedback (via SNR-induced shocks, 0,B star far-UV radiation) will only lenghten 
any such T* C h by destroying the H2 clouds fueling star formation (which then have 
to reform to be SF-capable again). From Fig. 11 . 61 it is obvious that for (U)LIRGs it 
will always be TAD>T CO ns, tnus rendering AD-regulated star formation as unrealistic. 
For the average ISM conditions in such vigorously star-forming galaxies a process 
such as reconnection diffusion ( ll36l ) can still provide magnetic field diffusion fast 
enough as to keep T <N . ch <T C ons (E2J). 

In summary, in the AD-diffusion regulated star-formation scenario, as in the 
gravoturbulent one, the large CR energy densities in the CRDRs of compact star- 
bursts will dramatically alter the initial conditions for star formation. The effects on 
the emergent stellar IMF are however less clear to discern than for gravoturbulently- 
regulated star formation, especially if fast reconnection diffusion plays a role in re- 
moving magnetic field support from dense H2 gas with orders of magnitude larger 
ionization fractions than in ordinary star-forming galaxies. 



1.4.2 CR-induced effects on the equation of state of the H2 gas 

Another effect of the large CR energy densities in the CRDRs of compact starbursts 
is the erasing of the so-called inflection point of the effective equation of state (EOS) 
of the H2 gas. The latter, parametrized as a poly trope P=Kp y , is used in many numer- 
ical simulations of self-gravitating turbulent gas. These found its polytropic index 
7, and more specifically the so-called EOS inflection point (i.e. the characteristic 
density n c where y< 1 forn<n c flips to y>\ for n>n c ), playing a crucial role in de- 
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1 2 3 4 5 

log [s- 1 ] 

Fig. 1.6 The ambipolar diffusion (AD) timescale as a function of the CR ionization rate 
C-17=Ccr/(10" 17 s- 1 ) (assumed ocUcr), for n=[10 4 , 10 s , 10 6 ]cm- 3 (n 4 =n(H 2 )/(10 4 cm- 3 )). For 
comparison the AD timescale for ISM conditions typical in the Galaxy is also shown as a horizon- 
tal line at the bottom. The tb gas consumption timescale T cons ~M(Ft2)/SFR for a typical (U)LIRG 
is also marked by the higher horizontal line. The shaded area marks the region where such galaxies 
are expected to be in terms of the £cr in their CRDRs, where clearly Tad > Icons- 



termining the functional shape and mass scale of the resulting mass spectrum of the 
collapsed dense gas cores (and thus the stellar IMF) ( ll38l l39l ). For ordinary ISM 
conditions this inflection point occurs within CRDRs at n c ^10 5 cm~ 3 ( Il22l [39l ) 
because of the onset of strong thermal gas-dust coupling in the (far-UV)/optically- 
shielded gas regions, which makes them nearly isothermal and sets a characteristic 
IMF mass scale of Mi ch ^Mj eans (n c ). 

In the CRDRs of compact starbursts the H2 gas and the dust remain thermally 
decoupled with Ti c ; n >T ( iust f° r densities as high as n^l0 6 cm~ 3 (Fig. I1.41 i. erasing 
the EOS inflection point. This maintains strong H2 gas cooling (^n 2 ) at much higher 
densities (>n c ), rendering local cooling timescales (t coo i) much shorter than the dy- 
namical ones (Tdyn)- This allows thermodynamics to always stay a step ahead of 
self-gravity in decreasing the local Jeans mass (as densities climb and temperatures 
drop) faster than the gravitationally-driven evolution of the cloud fragments it to 
mass fragments of ^Mj eans (r,t). A comparison of these two important timescales 
as a function of the CR energy densities is shown in Fig. 11.71 from where it be- 
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comes obvious that an inflection point where Tdyn/T CO oi^l * s not even approached 
once UcR/UcR,Gai>10, even for densities as high as ^10 6 cm~ 3 . Clearly any future 
numerical simulations of molecular clouds in extreme CRDRs must take this into 
account, or abandon the use of an EOS for the gas altogether. 



1.5 A CR-controlled stellar IMF: some consequences 

A top-heavy stellar IMF can be the result a higher characteristic mass M^. h for the 
young stars (the so-called IMF 'knee', which is ^(O.3-O.5)M0 for the Galaxy) or a 
shalower power law than its Salpeter value of a^2.35 (for an IMF functional form 
of dN(M*)/dM*ocM~ a and stellar masses between 1 M and 100 M Q ). In a gravo- 
turbulent scenario the exact shape and mass scale of the emergent stellar IMF under 
the very different star formation initial conditions in extreme CRDRs of ULIRGs 
can only be determined with new H2 cloud simulations that make use of these new 
conditions. In the past simulations using so-called starburst ISM conditions have 
been done, recovering large M^. h in such environments (|40|). Unfortunately such 
simulations used conditions typical of PDRs rather than CRDRs, an approach that 
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would have then produced a varying IMF even in the Milky Way given the large 
range of gas properties in its PDRs. It was only recently that the effects of large CR 
energy densities on H2 cloud fragmentation and the stellar IMF have been explored 
(ED) and found M^ boosted to -(2-3) M (Fig.[0]>. 




Fig. 1.8 The stellar initial mass functions (IMFs) from numerical simulations of turbulent H2 
clouds, for UcR=UcR,Gal (left), lOOxUcR.Gal (right) 1411 . For comparison puiposes, a Salpeter- 
type (green dashed) and Chabrier-type (blue dot-dashed) IMF fits to the emergent mass spectrum 
are also shown along with a linear fit and a lognormal fit (purple and red lines). The power-law 
slopes above the turn-over mass are given in the upper left corner. 



The dependance of the thermal and ionization state of the gas in CRDRs on the 
average Ucr permeating the ISM of a galaxy naturally links the SF initial conditions 
and the resulting IMF to the average SFR density p s f r , provided UcR^Psfr (with the 
details of CR transport/escape mechanisms in quiescent SF disks and starbursts set- 
ting the proportionality factor). CR-controlled star formation initial conditions then 
naturally yield a bimodal stellar IMF in galaxies, with high values of p s t r determin- 
ing the branching point, namely: 

(merger-driven star formation)— >(high p s fr)— Khigh Ucr)— ^(top heavy IMF), 
(isolated gas-rich disk star formation)— >(low p. v / r )— >(low Ucr)— >(Galactic IMF). 

Such an IMF bimodality, with a top-heavy IMF in merger/starbursts versus a regu- 
lar IMF in isolated gas-rich disks has been postulated for some time in hierarchical 
galaxy formation models seeking to explain their relative populations across cos- 
mic time in ACDM cosmologies ([42]). CR-controlled star formation initial con- 
ditions in galaxies can now set this on a firm physical basis for the first time. It 
must be noted that a top-heavy IMF in high-p 4 /, systems such as compact star- 
bursts in mergers, will also have serious implications on the SFRs deduced from 
their observed IR luminosities as such IMFs will have several times higher energy 
outputs per stellar mass than the ordinary one. Then the tremendeous SFRs deduced 
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for merger/starbursts, especially at high redshifts (^10 3 M Q yr _1 ), may actually be 
~(3-5) times lower. 

Finally an p^-dependant IMF as the outcome of CR-controlled star forma- 
tion initial conditions in galaxies will also naturally yield integrated galactic IMFs 
(IGIMFs) that depend on the star formation history (SFH) of galaxies. This is sim- 
ply because p s f r can change significantly during the evolutionary track of a galaxy, 
especially during its early gas-rich epoch, even in the absence of mergers. Such a 
SFH-dependance of the IGIMFs has been considered as the cause behind the well- 
known mass-metallicity relation of galaxies (|43|). 



1.6 Future outlook 

The unique signatures imparted on the thermal and chemical state of the H2 gas 
in extreme CRDRs have now been studied extensively ( ||29l l44l 145 1). and most of 
them will be easily accessible in the new era of ALMA whose commissioning is now 
ongoing at Llano de Chajnantor on the Atacama Desert Plateau. Numerical studies 
of H2 clouds without the use of an EOS for the gas, that take also into account 
the strongly anchored magnetic fields on dense gas with large ionization fractions 
are now becoming possible. High-energy y-ray observations now started directly 
probing the average CR energy densities of other galaxies ( fl46l [47]). Thus we are 
looking at a decade where the properties of CRDRs, their effect on the stellar IMF of 
galaxies, and their role in galaxy formation can be decisively explored with powerful 
theoretical and observational tools. Finally these explorations will critically benefit 
from studies of the stellar IMF and its dependence on ambient conditions in galaxies 
that are now converging to much sharper picture ([48 1). 
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Appendix 

The CR heating rate used in this work is given by 

r CR = 1.95 x KT 28 n H ( 13x j^i 7g -i ) ergscm-V 1 (1.10) 

where Ccr^Ucr being the Cosmic Ray ionization rate (in s ), with an adopted 
Galactic value of CcR.Gai=5x 10~ 17 s _1 (corresponding to UcR=UcR,GaiX an d n H=2n(H2) 
for fully molecular gas. For an optically thin OI line (and thus maximal cooling) 
Aoi63 = XonHQuE u i, which becomes 

Aoi 6 3 = 3.14xl0- 14 Zon H C u i^^e X p(-227.72/T kin )ergs (1.11) 

where g n =3 and g\=5, and Xo=[OfH] being the abundance of oxygen not locked onto 
CO (%o^4.89x 10~ 4 ). The collisional de-excitation coefficient is given by ( |fl6l ) 

Cui =n H 10 a321ogTk ">- 10 - 52 = 3.02 x 10- 11 n H T^ 2 cm- 3 s- 1 , (1.12) 
Thus finally Equation 1.11 becomes, 

Aoi63 =2.78 x 10- 28 n 2 1 TMn 2ex P(-227.72/T kln )ergscm- 3 s- 1 . (1.13) 

For densities of «h <10 5 cm -3 and strong CR fluxes, a small fraction of carbon 
remains in the form of C + , acting as a coolant with 

A c+ = 1.975 x 10" 23 n 2 1 ^ c+ exp(-92.2/T kin )ergscm- 3 s- 1 , (1.14) 

computed in a similar fashion as the OI(63 jim) line cooling, and for a fully trans- 
parent medium. Gas-grain accommodation cooling the gas depending on their tem- 
perature difference can be expressed as ( ifTTl ): 

A g _ d = 4.0x 10- 12 an H n d (7 gV /T~(T kin -T dust ) (1.15) 

where n g is the number density of dust grains and (7 g is the grain cross-section 
(cm~ 2 ). If the gas-to-dust ratio is 100, the dust mass density is Pd=3.5 gem -3 , and 
a dust radius of a = 0. 1 )J.m then 

n d = 2.2 x 1(T 2 x J Li H n H /(4/3^p d a 3 ) = 7.88 x 10~ 12 n H . (1.16) 
The gas-grain accommodation factor a is given by (| 18 1): 

a = 035cK V (- ] f^^\ (1.17) 
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which we set to the maximum value a=0.35 (i.e. maximum gas cooling from gas- 
dust interaction). Thus the cooling term due to the gas-dust interaction (Equation 
1.15) becomes 

A g -d = 3.47x lO-^v/Tk^^kin-Td^Oergscm-^- 1 . (1.18) 

Finally detailed PDR code models allow us to parametrize the cooling due to the 
CO rotational transitions as 

where we set £ C o/*[c]=(0.97,0.98,0.99,1.0) for densities of 5x 10 3 , 10 4 , 10 s , 10 6 
cirT 3 respectively (with Xco/X[c] = l corresponding to all carbon locked onto CO). 
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